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ABSTRACT 



We investigate the hypothesis that Lya absorption lines arise in two populations of 
halos — minihalos of small circular velocity (V^ ^ 55 km s~^) in which star formation 
On i and metal production are inhibited by photoionization of the UV background radiation, 

^— V I and large galactic halos (55 ^ Vc ^ 250 km s~^) which possess stars and metals. 

CN I Based on the model of Lya-absorbing gas confined in both populations of halos, we 

1^^ ' attempt to explain the recent observations of (1) associations of visible galaxies with 

On ■ Lya lines at low redshifts z ^ 1, and (2) metal lines associated with a non-negligible 

fraction of low H I column density Lya lines at z ~ 3. For galactic halos, we find 
that photoionized gas clouds confined in the pressure of ambient hot gas can produce 
2 . Lya absorptions with H I column density as low as 10^^ cm~^, and that the impact 

^ I parameter of a sightline for such absorptions matches well with the observed radius of 

gaseous envelope in a typical luminous galaxy. Using the Press-Schechter prescription 
for the mass function of halos, we also show that the fraction of Lya lines with 
^ ' associated metal lines can be understood in terms of the fraction of Lya absorbers that 

j^ ■ are associated with galactic halos. In particular, the reported fraction of ~ 0.5-0.75 at 

z ~ 3 is reproduced when the boundary value of Vc to separate mini or galactic halos is 
40 ~ 60 km s"'^, which is consistent with the theoretical prediction of galaxy formation 
under photoionization. The average metallicity of both Lya forest and damped Lya 
systems at z ~ 3 is explained in terms of the model of halo- formation history combined 
with the age-metallicity relationship of Galactic halo stars. Possible methods to test 
this hypothesis and the other alternative scenarios are also discussed. 



Subject headings: galaxies : abundances - galaxies : evolution - galaxies : intergalactic 
medium 
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1. INTRODUCTION 

Several recent observations indicate that even the low H I column density Lya absorption 
lines may not be primordial as were previously thought. High resolution observations by Cowie 
et al. (1995), Tytler et al. (1995) and Womble et al. (1996) show that there are metal lines (C IV) 
associated with Lya absorption lines with as small H I column density as A'^hi ^ 10"*^^ cm~^ at 
z ~ 3. The fraction of the Lya lines with N^j ^ 10^"^ cm~^ that are associated with metal lines is 
approximately ^ 0.5 at z ~ 3. The median abundance of lines with A^hi ~ 10^^'^ cm~^ at z ~ 3, 
is found to be [C/H] ~ —2.5. Recently Songaila and Cowie (1996) have shown that the fraction of 
Lya lines that are associated with metal lines, hereafter denoted as fm, could be higher, and they 
have also detected the other C II, S IV, N V absorption lines in hydrogen clouds with Nhi > W^^ 
cm^^. In view of the dependence of observational sensitivity for detecting metal lines, the fraction 
fm represents only a lower limit for the fraction of Lya lines associated with 'metals'. 

Another set of observations shows that at low redshifts a large fraction of Lya absorption 
lines are associated with galaxies (Maloney 1992; Morris et al. 1993). Lanzetta et al. (1995b, 
LBTW) claimed that ~ 60 it 19% of the Lya lines at z ^ 1 arise in galactic halos. The observation 
of Barcons et al. (1995) in which they imaged a galaxy corresponding to an absorption line, is an 
evidence of such associations. It is, therefore, reasonable to assume that at least a fraction of the 
Lya absorption lines arise in galactic halos, and this fraction is hereafter denoted as fg. 

In the hierarchical model of structure formation, there is a continuum of halos of various 
masses and velocity dispersions. In this paper, we consider the hypothesis that the Lya lines 
are produced in two populations of halos in this continuum - in 'minihalos' of small velocity 
dispersions which have not cooled and formed stars, and in large galactic halos with stars and 
metals. Based on this two-population model for QSO absorbers, we attempt to explain fm in 
terms of fg by the relation fm ~ fg- 

There are circumstantial evidences for such a two population scenario. In the column density 
distribution of Lya lines, there is a break around A'^hi ~ 10^^ cm~^ (see, e.g., Cristiani et al. 1995). 
Petitjean et al. (1993) have interpreted this break as being due to existence of two populations of 
Lya lines - the lines with higher H I column density being associated with strong metal lines and 
with galactic halos. Recently, Fernandez-Soto et al. (1996) analyzed the two-point correlation of 
very weak C IV absorption lines associated with high-redshift {z ~ 2.6) Lya absorption systems. 
They found that high-redshift Lya absorption systems traced by C IV lines are clustered in 
redshift, as strongly as that expected for galaxies. This suggests that many Lya absorbers even at 
high redshifts may arise in galaxies. 

Our primary goal in this paper is therefore to determine whether these recent observations 
can be understood within the framework of a two-population model for Lya absorption systems, 
including their redshift evolution in the context of hierarchical structure formation. In §2, we 
present a model of Lya-absorbing gas confined in both mini and galactic halos and set the 
prescription for determining the number of absorption lines provided by various masses of halos. 



On the basis of the Press-Schechter mass function with various power spectra and cosmological 
parameters, we then discuss the fraction of absorption hues in minihalos and in galactic halos by 
comparing with the recent Keck observations in §3. The simphfied approach for modehng chemical 
evolution of halos is also presented to predict the evolution of the average metallicity associated 
with Lya lines in redshift. The discussion and conclusions are finally drawn in §4. 



2. MINIHALOS AND GALACTIC HALOS 

2.1. Assumptions and definitions 

We define minihalos as halos in which baryonic gas is entirely photoionized by the external 
UV background and the gas is stably confined in the gravitational field of the dark matter (Rees 
1986; Ikeuchi 1986). In these minihalos, the gas does not go through extensive star formation 
because of suppression of radiative cooling by photoionization (Efstathiou 1992). Below, we will 
characterise the halos by their circular velocities Vc, where the mass density of halos is dominated 
by extended dark matter with r^^ distribution. Minihalos are therefore bounded by two scales of 
circular velocity, Vi < Vc < V2. The lower limit is set by the sound velocity of the photoionized 
intergalactic medium (IGM), which is ~ 15 (T/gm/IO^ K)"'^ km s~^. Halos with Vc less than this 
value will not have baryonic gas infall from the IGM. This is basically the Jeans mass, below 
which the perturbation in baryonic gas is suppressed. We take a characteristic value of Vi = 15 
km s~^ for the photoionized IGM with Tjgm ~ 10*^ K- 

The upper limit V2 is admittedly less certain. Recent numerical simulations show that 
photoionization severely inhibits the formation of galaxies with Vc ^ 30 km s^^, and substantially 
decreases the mass of cooled baryons in halos with Vc ^ 50 km s^^ (Thoul and Weinberg 1996). 
They also find that above Vc ~ 75 km s~^ photoionization of the IGM does not have any significant 
effect. For the reasons described in the next subsection, we will use a canonical value of V2 = 55 
km s~^ and examine the effect of changing V2 over its possible range. In fact, we later ask the 
question whether or not one can determine the value of V2 from the observations. In brief, then, 
we assume that halos with Vc < ^ ~ 55 km s^-*^ do not undergo vigorous star formation, and call 
them minihalos in this paper. 

We will call the halos with circular velocity larger than V2, the galactic halos. Our results are 
not sensitive to the upper limit of the circular velocity for galactic halos, and we will use a value 
of V3 = 250 km s~^ as an upper limit. Below, we will associate these halos with Lya lines which 
have corresponding metal lines. It is true that the division of halos into the above categories of 
mini and galactic halos is rather ad hoc. However, it is most likely to be close to the real picture 
for the above mentioned reasons. 

For the gas associated in galactic halos, the virial temperature T for I4 ^ V2 ~ 55 km s^^ is 
above 10^ K [T ~ W^{Vc/55kms~^)'^ K], where the mean molecular weight is taken as 0.59 for 



a totally ionized gas. In this temperature range, any local disturbances in the hot gas may be 
thermally unstable, because the radiative cooling rate increases with decreasing T. We therefore 
assume that the instability promotes the development of a two-phase structure in the gas (Fall &: 
Rees 1985), so that it is the cool gas clouds, being pressure-confined by the surrounding hot gas, 
that produce Lya absorptions, as proposed by Mo (1994). In this work, we will be concerned with 
Lya lines with low HI column density as A^hi ~ 10^^ cm^^, in which the shielding against the UV 
background radiation is ineffective. The internal density of clouds in the outer parts of halos that 
give rise to such lines is low enough to be photoionized by the UV background radiation. We thus 
assume the gas in clouds to be in photoionization equilibrium, where the gaseous temperature 
is expected to be order of few x 10^ K, and therefore consistent with the observed Doppler 
parameters of Lya lines (see, e.g., Maloney 1992; Ranch et al. 1996). 

Note that, even the gas in disks, if disks are formed, will be ionized by the UV radiation in 
the outer parts, beyond the observed sharp edge at A'"hi ~ few x 10^^ cm^^ (Maloney 1992). The 
galaxies at low redshift imaged by Barcons et al. (1995) are inferred to give rise to absorption 
lines at an impact parameters of ^ 50/i^^ kpc, where Hq = lOO/i km s~^ Mpc^^. The gas which 
is responsible for such absorption lines is certainly ionized by the UV background radiation (see 
also, Morris et al. 1993; Salpeter and Hoffmann 1995). 



2.2. H I column density 

We assume that the density profile of both mini and galactic halos is represented by that of a 
softened isothermal sphere, 

P=^^^ (1) 

j,Z _j_ j,Z 

where p is the central density and Tc is the core radius. The baryonic gas amounts to a fraction F 
of the total mass of the halo. Then, the circular velocity Vc{r) at r can be written as. 



GM(r) 1 

V^{r) = — = ATTGprl-{x - arctanx) , (2) 

r X 



where x = r/r^- This functional form of V^(r) assures a nearly flat rotation curve at large radii 
r S> Vc- In what follows, the circular velocity is characterised by its value at the virial radius of 
the halo, r^j^ (see below for an explanation of r^j^)- Then, substituting eq.(2) into eq.(l) and 
introducing x^ = rmr/fc^ the density profile of a halo that includes both baryonic gas and dark 
matter is re- written as, 

4:-KG{r^ + r1) ^Xv — arctanxt, 
where Vc is the circular velocity at r = r^j^- 



P(r,v^) = ^z?T7±—rA—^^:rzrz-) ^ (3) 



We model the gas in minihalos {Vi < T4 < V2) to be in photoionized equilibrium with the 
external UV background radiation. We write the intensity of the UV background radiation at the 



Lyman limit as J = J_2il0^^^ erg cm^^ s^^ sr^^ Hz^^. Using the ionization and recombination 
rates for a fiducial UV spectrum of Jy oc u^"^ (Black 1981), the H I number density uhi in 
minihalos is written as, 

^Hl = -T—j ^n-2lr< \ J ^°^ mmihalos , (4 

where Gh = 2.54 x 10^ for hydrogen and an is the recombination rate coefficient (Black 1981). 
The column density of H I at an impact parameter b is then estimated as 



""'<" = 'I 7P^* 



14/ T \-3/4. F \2 _-^ / Vc 



cm for minihalos . (5) 



Xj, — arctanXj,/ [(6/10fcpc)2 + (rc/lOfcpc)^]^ 

For galactic halos (V2 < K: < V3), we assume a cloud confined by ambient hot gas to be a 
uniform sphere with proton number density rici and temperature T^i- The pressure balance with 
hot gas with proton number density n (= Fp/nip) and temperature T yields rid = nT /Td- Under 
the assumption that the clouds are photoionized by the external UV background and in ionization 
equilibrium, the H I number density nui is then written as 

UHI = -. — J -I t-^t^ ^ J-2i'^H{Td)(7^) r? for galactic halos . (6) 

The cloud radius Rd for the mass of each cloud Md is 

yAnndrnp^ " \Td'' 

Each cloud produces the H I column density as NHi{r) ~ nHiir)Rdi'>')- In order to obtain the H I 
column density observed at an impact parameter b from the galactic center, it is required to know 
the spatial distribution of clouds in halos a priori. Without going into details for modeling this, 
we assume that the clouds have roughly unit covering factor in halos when projected onto the sky, 
and that each sightline encounters only one cloud in average, which is not in disagreement with the 
observational results (Morris et al. 1993; LBTW). Then, as we argue in the Appendix, the average 
H I column density observed at an impact parameter b is well approximated by Nniir -^ b). We 
thus obtain from eq.(6) and (7) 

AT ^A^ n« -,nl5{ ^cl yrof Td ^-29/12 f F ^5/Z f V^ ^20/3 

5/3 1 _2 



arctan Xy > [(VlOOfcpc)^ + (rc/100A;pc)2]5/3 



cm for galactic halos ,(8) 



where we assume the virial temperature for T of the hot gas. We note that both expressions for 
^Hl given in eq.(5) and (8) are approximately connected continuously at the canonical value 



of 14 ~ 55 km s^^ for V2, corresponding to the virial temperature of T ~ 10^ K; this is the 
temperature of hot gas below which the picture of a two-phase structure in the gas becomes 
invahd. We therefore use eq.(5) for Vc < 55 km s~^ and eq.(8) for 14 > 55 km s~^ in the following. 

According to several determinations of the intensity of the UV background radiation (e.g. 
Madau 1992), J-21 is nearly constant at high redshifts 2.5 ^ z ^ 3.5, whereas at low redshifts 
z ^ 2.5 it rapidly decreases with decreasing redshift. We adopt the following functional form for 

J-21, 

J_2i for z > 2.5 



'-''^'^ = { j-4^y for.; 2.5, ^'^ 

where the index a quantifies the decreasing rate of J_2i for z < 2.5. 

To estimate the value of Xv{= Vmr/rc) in eq.(3), we use the theory of the spherical collapse 
of cosmological density perturbations, that yields r^r = O-^Hq (1 + z)^^"T4 for the Einstein-de 
Sitter (EdS) universe (e.g. Peebles 1980; Padmanabhan 1993). This allows us to derive x„ using 
the observed values of Vc and Tc of the present-day dark matter halos at z = 0. For this purpose, 
it is convenient to use dark-matter dominated dwarf spirals (Moore 1994), where the dynamical 
effect of a luminous component on the observed rotation curve is sufficiently small, so that one 
can accurately determine the structural parameters of dark halos. Adopting the values of Vc and 
re from Moore (1994) for dwarf spirals of DD0154, DDO170, DDO105 and NGC3109, we obtain 
W ^ x^ ^ 20 for the EdS universe with h = 0.5. In what follows, we assume x^, = 15 = const., 
since the following results in the concerned low-column density range of Nhi rS 10^^ cm~^ are 
quite insensitive to the values of x^ as long as x^, ^ 1. 

Figure 1(a) shows the impact parameter b vs. the total mass of halos M at z = 0.5, where we 
take r = 3 X 10"^ K for the gas in minihalos (eq.5), T^ = 3 x 10*^ K and Md = 10^ Mq for the clouds 
in galactic halos (eq.8), with the common parameters of F = 0.05, J-21 = I, a = 2, and x^ = 15. 
It is seen that lines of constant Nhi have the slightly different slope at M ^ 2.3 x 10^*^ Mq (or 
Vc ^ 55 km s""*^) from that at larger M or Vc. This is because Nhi oc b~^ in eq.(5) while oc b~^^'^ 
in eq.(8) for b ^ Vc. Here we note that the discontinuities at the boundary of both regimes (at 
T4 ~ 55 km s^^) are quite small in b so that these will not affect the later results. We also find that 
H I column densities as low as Nhi ~ 10^^ cm^^ can be provided by either the low-mass minihalos 
with the impact parameter b below ~ 40 kpc or by the high-mass galactic halos with larger b. 
This is more clearly presented in Fig. 1(b) that shows Nhi vs. b for a given total mass M. Dotted 
line denote the observed size of gaseous envelopes ~ 160 h^^ kpc in luminous galaxies which give 
rise to Lya absorptions at z ^ 1 (LBTW). In our model of absorbers, this is achieved for a mass 
range of 10^^ ^ M ^ 10^^ Mq, which corresponds to galaxies with luminosity 10^ ^ L ^ 10^*^ Lq 
for the baryon fraction F of 0.05 and the mass-to-light ratio of the order of 5. Thus, the model 
yields the right order of the size of extended Lya-absorbing gaseous envelopes in galaxies. 

Bold solid line in Fig. 1(b) corresponds to the observed correlation between the H I column 
density and impact parameter by LBTW, assuming Doppler parameters of absorption lines in the 



range of 20 — 40 km s^^. The present model also implies the correlation and agrees approximately 
with LBTW, where the slight discrepancy in the slope may be due to the effects of sampling over 
different galactic masses and/or impact parameters. Recently Le Brun et al. (1996) and Bowen 
et al. (1996) claimed that the equivalent width of Lya lines do not correlate with the impact 
parameter. This has been argued to weaken the possibility of Lya lines being directly associated 
with galaxies. The other possibilities discussed in the literature include the scenario of filaments 
between luminous galaxies, and Lya lines from outflows from galaxies. The issuse may be settled 
by assembling the statistically meaningful number of Lya absorber-galaxy pairs and by comparing 
with the current type model that presents the explicit dependence of an impact parameter on the 
masses of galaxies and redshifts (Fig.lc). 



2.3. Cosmological context 

With the knowledge of the H I column density that is produced at an impact parameter 
b{Niii, Vc), we can calculate the number of Lya lines with a column density larger than N^i, if the 
number density of halos with circular velocity Vc is known. We use the Press-Schechter formalism 
for this purpose, and CDM models of structure formation (similar to the work of Mo et al. 1993). 

For the EdS universe, the circular velocity Vc and mass M of a CDM halo are related as 

A'jr 

M = —porl, Vc = l.67{l + zf-^Horo, (10) 

where po is the present-day mean density of the universe and vq is the comoving radius of the halo. 
The comoving number density of halos per unit circular velocity Vc is given by (Mo et al. 1993) 

- -3(L67)3j.go(l + .)V^dlnA 
n[Vc,z)aVc - (2vr)3/2iMA(ro) dlnVc 

Here A(ro) and 6c = 1.68 are the rms and threshold linear overdensities in a spherical region of 
radius ro, respectively, extrapolated to the present epoch. With a top-hat filtering, the functional 
form of A(ro) for the CDM power spectrum of density perturbations is taken from Bardeen et al. 
(1986), and its value is normalized at a scale of 8h^^ Mpc with the bias parameter bg. 

The number of absorption systems per unit redshift with an H I column density larger than 
A^'hi, that is produced in halos with circular velocity Vi < Vc < Vu, is then given as (Sargent et al. 
1980; Mo et al. 1993), 

dNiN^„z,Vi,Vu) c 1/2 r^ ,,r t^r \ h^(AT T/^ (^^^ 

= —-(1 + 2;)/ / dVcen{Vc, z)T:b {N^^,Vc) , (12) 

dz Hq Jvi 

for the EdS universe. We define {Vi,Vu) = (Vi,V2) for minihalos and (V2,V3) for galactic halos. 
Here e denotes the fraction of halos that give rise to Lya absorptions, in view of the fact that 



some fraction of halos, specifically early-type galaxies such as E, SO, dE and dSO in the scale of 
galactic halos, are lacking of interstellar gas. Since these galaxies experienced the burst of star 
formation and the subsequent gas loss via galactic wind for only a short period of rS 1 Gyr (see 
e.g. Yoshii & Arimoto 1987), the probability of a sightline to pass the protogalactic, gas-rich stage 
of these galaxies may be small. To evaluate e, we adopt the nominal fraction 69% of late-type 
galaxies among all types derived by Postman and Geller (1984) and assume it to be constant with 
z for simplicity. Thus, 

{1 for mini halos 

0.69 for galactic halos . 

Then the above equation can be used to find the fraction of Lya lines with > A'^hi that are 
associated with galactic halos: 

^ galactic _ d]SS{N^^,z,V2,Vz)ldz 
^ mini -I- galactic dN{NYii^z^V\^V-i)ldz' 

Besides the case of the 'standard' EdS universe, we will also calculate the number of Lya 
lines dN/dz and their galactic fraction fg for a world model with a cosmological constant. The 
basic cosmological parameters that define the evolution of the universe are the density parameter 
i7o = Pq/ Pci where pc is the critical density to close the universe, and the cosmological constant 
Aq = Ac^ /3Hq. We investigate a low-density (JIq < 1); flat {Qq -|- Aq = 1) universe, since such a 
universe deserves special attention to solve several paradoxical results of cosmological observations 
(Ostriker & Steinhardt 1995). We use the corresponding set of equations (10)-(12) derived from 
the spherical collapse theory of density perturbations in a non-zero Aq universe (see Peebles 1980 
and Suto 1993 for details). The model universes which we adopt for the following analysis are 
{Qq, Aq, h) = (1, 0, 0.5) and (0.4, 0.6, 0.65). Both universes yield the cosmological age as ~ 13 
Gyr, so that any differences in the final results, which may arise from the different cosmological 
age, are eliminated. 



3. MODEL RESULTS 

3.1. The number of absorption lines 

We first calculate the number of Lya absorption lines based on the model described in the 
last section, and investigate whether or not the model results can be reconciled with observations. 

Figure 2(a) shows the redshift evolution of the number of Lya absorption lines with an H I 
column density larger than Nhj = 10^^ cm^^. The absorption lines arise from both minihalos and 
galactic halos, and the bounds on circular velocities are taken as Vi = 15 km s~^ and V3 = 250 
km s~^. We take a standard set of parameters for photoionized gas (T = 3 x 10^ K for minihalos, 
Tc/ = 3 X 10^ K and M^ = W^Mq for galactic halos, with F = 0.05 and J_2i = 1). Sohd fines 
correspond to a 'standard' EdS universe {Qq = 1, Aq = 0), while dotted lines to a low-density 



(r^o = 0.4), flat (r^o + Ao = 1) universe. It is found, as was already argued by Mo et al. (1993), 
that there are two stages in the evolution of the number density. First, the rapid formation of 
minihalos inherent in the CDM models gives rise to the abundant of Lya lines at z ;^ 10 — 20, and 
second, at lower redshifts of z ^ 10, the number of Lya lines decreases with time as the hierarchical 
mergings proceed — small-sized, less-massive objects disappear by being incorporated into larger, 
more massive objects. The rapid decrease of dN/dz with decreasing z is now an observationally 
established fact, first recognized by Peterson (1978). For the cases of large hg and non-zero Aq, 
the maximum value of dN/dz is reduced and the epoch at which it is realized is delayed. This is 
because the growth of density perturbations and subsequent mergings are suppressed due to the 
bias in galaxy formation and the rapid expansion of the universe. However, the effects of these 
different assumptions on model parameters appear to be almost indistinguishable at the second 
stage of the number evolution of Lya lines. This indicates that probing world models by dN/dz 
of absorption systems in the observed redshift region of z ^ 3 may not be useful in view of some 
model uncertainties. We note here that numerical simulations with a variety of world models have 
shown almost similar H I column density distribution for Lya lines (Miralda-Escude et al. 1995; 
Hernquist et al. 1995). Also plotted in the figure as filled squares are the observed numbers of 
Lya lines at z ~ 3 (Petijean et al. 1993), z ~ 1.5 (Lu et al. 1991), and z ~ (Bahcall et al. 1993). 
It is clear that the models with J-21 ~ 1 = const, at z ^ 2.5 are reasonably in agreement with the 
observed values in the corresponding range of redshifts, whereas the Hubble Space Telescope (HST) 
result of the unexpectedly large number of low-redshift Lya lines (Bahcall et al. 1993), compared 
to the extrapolated value from high z, can be explained if the UV flux is decreasing with time (Mo 
et al. 1993). We find that J-21 = [(1 + z)/3.5]° with a ~ 2 fits well. It is interesting to note that 
with such a redshift evolution, the current UV flux is J-2i{z = 0) ~ 8 x 10^^, which is consistent 
with the observations of Maloney (1992) and Kulkarni and Fall (1993). 

The dependence of dN/dz on the column density Nhi is shown in Fig. 2(b), where filled 
squares are taken from Ranch et al. (1992) for 2 ^ z ;$ 3. The observed power-law decline of 

—5/3 

the Lya lines, oc Njjj , is well reproduced by the model on the basis of the isothermal density 
profile of minihalos (Rees 1986), together with the current model of pressure-confined clouds in 
galactic halos. We remark that this is the case as long as the effect of self-shielding against the 
external UV radiation is neglected. The effect, leaving a cooled neutral core in the inner part of 
a sphere (e.g. Chiba Sz Nath 1994), is realized at the high column-density ends of Nhi ~ 10^^^^^ 
cm^^. This effect is thus unimportant in the present study of low column-density Lya lines of 
10^^ ^ Nhi ~ 10"*^^ cm^^ that occur at the large impact parameters of sightlines (Fig.l). 



3.2. Fraction of Lya lines associated with galactic halos 

We now investigate what fraction of Lya lines can be associated with galactic halos where 
production of metals is possible — the fraction fg which may be compared to that of Lya lines 
associated with metal lines, fm- 
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Figure 3(a) shows the evolution of the galactic fraction fg of Lya lines with Nhi > 10^^ 
cin~^, where we assume V2 = 55 km s^^ as the boundary value of Vc between minihalos and 
galactic halos. It is found that initially at high redshifts (z ^ 20), most of Lya lines are provided 
by minihalos that collapsed early for the CDM-type power spectrum of density perturbations. 
This leads to the small values of fg less than 0.2. Then, at 10 ^ z ^ 20, the fraction of lines in 
galactic halos is rapidly increasing with time as a result of the formation of massive halos via 
hierarchical merging. Finally at z ^ 10, the value of fg slowly converges to the current fraction 
of 0.75 ~ 0.80. The figure also indicates that the differences arising from the different bias 
parameters bg and cosmological parameters (rio,Ao), which are clearly visible at high redshifts, 
become less significant at lower redshifts where observations of absorption lines are accessible. All 
of the predicted galactic fractions /^ at 2; ~ 2.5 are in between the reported values of fm for metal 
lines, fm ~ 0.5 - 0.6 (Cowie et al. 1995; Tytler et al. 1995; Womble et al. 1996) and fm ~ 0.75 
(Songaila and Cowie 1996). 

The dependence of fg on the column density Njji at z = 2.5 is shown in Fig. 3(b). There is an 
indication of the increase of fg with Nhi', it would tend to be 1 at Nhi around lO^*^^^^ cm^^ for 
damped Lya clouds. There are no observational support yet for the above behaviour of fg with 
Nh[ and will be a test for the two population model. Recently, Songaila and Cowie (1996) claim 
that fm ~ 0.9 for Nhi > 1.6 x 10^^ cm'^ and fm ~ 0.75 for Nhi > 3.0 x 10^"^ cm'^, whereas 
the other observations give a much lower value of fm- More observations are needed to settle this 
issue (see also discussion). 

As we mentioned in §2, the exact value of V2, which separates the range of circular velocities 
into those of minihalos and galactic halos, is yet to be settled from theoretical points of view. 
For example, even for Vc ^ 55 km s~^ where the virial temperature is T ^ 10^ K, some fraction 
of gas in an inner part of minihalos may be able to cool and form stars, because of self-shielding 
effects against the external UV background radiation. We thus turn the argument and address the 
question, as to which range of V2 is allowed in the light of the Keck results. Figure 4 shows the 
galactic fraction of Lya lines as a function of V2, while Vi = 15 km s~^ and V3 = 250 km s~^. Here 
we note that the change of V2 is applied to the range of integration in eq.(12) and the definition 
of e in eq.(13), whereas the range of Vc for the expressions of Nhi (^ < 55 km s^^ for eq.5 and 
Vc > 55 km s~^ for eq.8) is unchanged on the theoretical grounds. For the observed fraction of 
0.5 ~ 0.75, the possible value of V2 to reproduce it is 40 to 60 km s~^. This is rather insensitive 
to the specific assumptions on the bias parameter bg, cosmological parameters (rio,Ao), and the 
limiting circular velocities (VijV^). However, to set a more robust constraint on the range of 
circular velocities or masses for mini/galactic halos, further observational determinations of the 
fraction of Lya lines associated with metal lines are needed. 
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3.3. Metallicity of absorption systems 

Several recent observations have also measured the metal abundances of absorption lines 
at high redshifts (Cowie et al. 1995; Womble et al. 1996). The presence of heavy elements 
implies that interstellar gas confined in dark halos has been processed by star formation and 
subsequent chemical evolution. While there are some alternative possibilities for enriching gas 
— one representative idea is that IGM as a whole has already been processed on a sub-galactic 
scale prior to formation of collapsed galactic halos (Songaila & Cowie 1996), we present here the 
hypothesis that enrichment occurs after halo formation. We postulate that gas within halos has 
been cooled quickly enough to fragment into stars, and being enriched. In our model, this is 
applied to the gas in galactic-sized halos with V2 < Vc < Vs, where cooling time is well less than 
dynamical time (Dekel & Silk 1986; Efstathiou 1992). Here, we attempt to estimate the possible 
metal abundances of gas in these halos. 

Suppose that a halo with Vc at the redshift z is enriched with an average metallicity ZiVc, z). 
Here the dependences on Vc and z arise from the fact that the halos with different Vc collapse and 
form at different epochs, say z/, leading to different epochs of star formation, say zsf^ and that 
the subsequent chemical evolution starting at zsf yields the metals at the redshift z. Thus an 
ensemble of these halos with the metallicity ZiVc, z) contribute to the metal abundance of Lya 
lines. Then since the probability for a sightline to encounter galactic halos in the velocity range of 
V2 <Vc <Vs IS given by eq.(12), we arrive at the average metallicity of Lya lines with H I column 
density larger than Nhi: 

< Z > (> Nhi, z) = -^ — T7 . (15) 

f^^dVc-n{Vc,z)T:b^{N^,,Vc) 

The actual expression for Z{Vc, z) in each halo is admittedly not straightforward to derive 
and very model-dependent. Physics involved in it are, e.g. merging history of a population of 
dark halos, gasdynamical and thermal processes of gas confined in these halos, and history of 
star formation. While the extensive and sophisticated approaches of incorporating these processes 
exist (Lacey &: Cole 1993; Kauffmann et al. 1993), we adopt a more simplified approach which 
we believe includes the essential properties of the system to deduce the characteristic values of 
metal abundance. First, we assume that the formation redshift Zf of the halo with mass M 
is represented by the collapse epoch of a rms, 1 a linear density contrast, i.e., if D{z) denotes 
the growth rate of density perturbation, Zf is derived from D{0)/D{zf) = A{M)/6c, where M 
is related to Vc from eq.(lO). According to Lacey & Cole (1993) on the basis of their extensive 
analysis of merging halos, this scaling relation appears to hold for the typical formation epoch 
of the main parent halo which had half or more of the present mass M (see their Fig. 9 and 
Fig. 10). Realistically the halo formation epoch for a given M is distributed with some probability 
arising from different amplitudes of density perturbations, but it can be characterised by the 
collapse epoch of a typical, 1 a perturbation. Second, gas in a virialized, galactic-sized halo 
must cool quickly (within a dynamical time) to ensure star formation, so that zsf, the epoch 
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of star formation, is assumed to be equal to Zf. This is reasonable, because the cooling time 
in the galactic scales is less than 10^ yr, which is a negligible time span — much less than the 
Hubble time — in the cosmological context. Third, concerning the subsequent chemical evolution, 
we adopt the observed age-metallicity relation of metal-poor halo stars in our Galaxy (Schuster 
& Nissen 1989), on the grounds that typical sightlines to produce the Lya lines with low H I 
column density penetrate the halo gas, from which old halo stars may have formed. The average 
enrichment rate AZ/At derived by Schuster & Nissen (1989) is of the order of 10^^ Gyr^-*^ for 
their sample of halo stars. This value is also close to the value from the chemical-evolution model 
of the halo by Pagel (1989). Therefore, denoting the age of the universe at the redshift z is tage{z), 
the metallicity Z at z is derived as, 

rtageiz) /^Z 

ZiV,,z)= —dt, (16) 

Jtage{zsF{Vc,z)) At 

where following the above arguments, the epoch of star formation, tsF ^ iage(-25F(V'c, z)), depends 
on Vc and z. 

We assume here that at the outer parts of the halo, which concerns us here for the Lya forest 
lines, the gradient of metallcitiy is negligible. This is a reasonable assumption because there are 
indications that star formation in the outer halo of our Galaxy has been slow and, therefore, the 
chemical evolution (Matteucci and Frangois 1992). If the variation in star formation rate is caused 
by the variation in the gas collapse time, then this variation is expected to be small in the outer 
regions of the halo. This assumption is consistent with the observation that the average metallicity 
is almost constant for Lya lines with Nhj ^ 10^^ cm~^ (Cowie et al. 1995, Songaila and Cowie 
1996). 

Figure 5(a) shows the evolution of the mean abundance for absorption lines with Nhj > 10^^ 
cm~^ (solid and dotted lines). It is found that for 3 ^ z ^ 5, the metal abundance significantly 
increases by more than two orders of magnitude. At z ~ 2.5, it is approximately in agreement 
with the Keck results of [C/II]= —2.5 ~ —2. We remark that the adopted age-metallicity relation 
of a constant enrichment rate overestimates the halo abundance at z ^ 1: the relation should be 
applied to stars with metallicity less than 10^^ of solar, which is the criterion to separate halos 
stars from disk stars (Schuster & Nissen 1989). The subsequent formation of halo stars is stopped, 
and most of the enriched halo gas would fall into the disk — the hypothesis supported for solving 
the so-called G-dwarf problem in the solar neighborhood (e.g. Yoshii et al. 1996). Thus we assume 
log < Z > /Zq = —1 as an upper limit for abundance of absorption lines in Fig. 5(a), but this 
does not change our result in the light of the Keck results at z = 2 ~ 3. We have also investigated 
the case when the epoch of star formation tsF is constant, being independent of the mass scale 
of halos — as was done by Timmes et al. (1995) for damped Lya systems. We find that the 
evolution of < Z > based on our prescription of star formation beginning at the collapse epoch of 
1 a density contrasts is essentially the same as the case when tsF = 1.5 ~ 2 Gyr. 

Also shown in Fig. 5(a) with filled circles are the recent findings of Lu, Sargent &: Barlow 
(1996) for the abundances ([Fe/H]) of damped Lya systems observed with the Keck. Intriguingly, 
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the median values of abundances in such systems are well reproduced by the current models. This 
may support the hypothesis, first proposed by Lanzetta, Wolfe &: Turnshek (1995a), that damped 
Lya systems could represent a spheroidal component of galaxies (see discussion in §4.2). 

The dependence of < Z > on the value of V2 is shown in Fig. 5(b). In contrast with the 
galactic fraction fg of Lya lines, < Z > turns out to be sensitive to both the bias parameter bg 
and the cosmological parameters (rio,Ao) at large values of V2. It is systematically lower for a 
larger bias and/or lower-density universe. This is caused by delayed star- formation epoch tsF, 
leading to less time for enriching gas. 



4. DISCUSSION AND CONCLUDING REMARKS 

4.1. Evolution of galactic fraction and average metallicity 

We have shown that the observed fraction of Lya absorption lines associated with 
metal lines can be understood in terms of photoionized gas confined in galactic-sized halos with 
40 — 60 ^ T4 ;$ 250 km s^^, whose number density follows the Press-Schechter model of hierarchical 
structure formation. The derived average abundance of heavy elements in absorption systems, 
combined with the characteristic epochs of star formation and the observed age-metallicity relation 
in our Galactic halo, appears to be reasonably in agreement with the observed abundance at 
z = 2 ~ 3. 

The current model also suggests that both the fraction of Lya lines associated with galactic 
halos and the metallicity in absorption systems are expected to increase with decreasing redshift 
(Fig. 3a and Fig. 5a). While the explicit dependences of these quantities on redshift are yet to 
be settled (see Songaila &: Cowie 1996), we note the findings of Steidel (1990) that the number 
density of C IV lines increases with decreasing redshift. In our model, this can be explained from 
the above predicted properties of absorption systems in the following way. 

In photoionization equilibrium, the ratio between the column densities of C IV and H I is a 
function of the ionization parameter F, which is the ratio of the densities of ionizing photons and 
particles, and of the metallicity Z; using the photoionization code of CLOUDY (Ferland 1993), 
we find that for log^^Q F = —1.5 to —2.0, and for Z = 0.001 to 0.5, the ratio Nciv/^Hi is well 
described as, 

^^-"'^^_1.3±o.2^1.0±0.1_ (17) 



Nhi 

Then, suppose that the number density of Lya lines per unit H I column density and redshift, 
d^Nni/dNnidz, is represented by a power law form. 



d'Nm 
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dNHidz(xN^i{l + zfdNHidz, (18) 



dNnidz 
where (5 ~ 1.5 (Sargent et al. 1989) and 7 ~ 2.45 (Lu et al. 1991), and that a fraction fg of these 
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Lya lines are associated with C IV lines, which may be written as, 

fg{z)^{l+Zr^. (19) 

From Fig. 3(a), we obtain ij ~ 1.6 at z = 2 ~ 3. 

Suppose also that the mean metallicity of absorption lines increases in time following 
< Z >oc (1 + z)~ . Then eq.(17) and eq.(18) give the column density distribution of C IV lines as, 

"^ ^^^^ dNcivdz oc iV^fy(l + zV^dNcivdz , (20) 



dNcivdz 
where 



7C = 7 - S{P) - 11 , (21) 



describes the redshift evolution of C IV lines. It is known from observations that the [3 index is 
the same for both H I and C IV lines (Cowie et al. 1995; Bergeron et al. 1994), indicating that 
the simple equation (17) holds well. 

Steidel (1990) has shown that the index 7c is negative, of order —1.25. This means that the 
number of C IV lines increases in time. In our model, this can be due partly to the increase of the 
galactic fraction fg in time and to the increase of the metallicity < Z > in time. Then, for the set 
of power-law indices, f3 = 1.5, 7c = —1.25, 7 = 2.45, and ij = 1.6, we arrive at 6 = 1.4, thereby 
indicating, 

< Z >oc {1 + zr'^-\ (22) 

It is worth noting that the enrichment rate of absorption lines shown in Fig. 5 (a) is well 
approximated by this equation. The redshift evolution for C IV lines used here is, strictly 
speaking, for high C IV column density lines. It will be interesting to see whether observations in 
the near future show a similar redshift evolution for small C IV column density lines. 

Therefore, the derived redshift evolutions of both the fraction of Lya lines associated with 
galactic halos (Fig. 3a) and the metallicity in absorption systems (Fig. 5a) are approximately 
consistent with the observed scaling laws for Lya and C IV lines. 



4.2. Possible relation with metals in damped Lya systems 

The abundances of metals have also been measured in damped Lya systems, which are 
believed to be protogalaxies. Smith et al. (1996) and Pettini et al. (1995) have recently reported, 
by measuring the Zn lines as a probe of metallicity, that the typical Zn abundance at z ~ 2 is 
1/15 of solar, but there is a considerable scatter in the Zn abundance, spanning more than two 
orders of magnitude. Using high quality spectra obtained with the Keck telescope, Lu, Sargent 
& Barlow (1996) have also derived a large scatter of abundances, about a factor of 30, with 
-^///-weighted mean abundance of 0.028 solar between 2 < z < 3. In the context of the present 
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model, such a large scatter may be caused by the different formation epoch of halos, Zf, for a 
given mass arising from the different amplitudes of primordial density fluctuations and different 
merging history (Lacey & Cole 1993). Also, for the adopted age-metallicity relationship of halo 
stars, there is a scatter in the ages of the stars of ±2.5 Gyr at a given abundance, leading to a 
scatter of abundances more than an order of magnitude at a given age (Schuster & Nissen 1989). 
It is worth noting that the evolution of the mean abundance is reproduced from the formation 
epoch of a typical 1 a fluctuation and the mean age-metallicity relationship of halo stars in our 
Galaxy (Fig. 5a). This suggests that damped Lya systems represent a metal-poor halo component. 
Also, Songaila and Cowie (1996) found that the ratio Si/C is about three times solar, similar to 
the abundance ratio of Galactic halo stars. 



4.3. The possibility of an enriched IGM 

The two population model, as presented here, clearly predicts (a) the fraction of Lya lines 
with associated metal lines as a function of redshift and (b) the redshift evolution of the average 
metallicity. The observations, however, have not yet converged to any firm conclusion and it is 
too early to rule out other possibilities. As noted in Songaila and Cowie (1996), it will be very 
difficult to determine the fraction of metal enriched Lya clouds below Nhi ~ 10^^ cm~^, because 
this will mean a sensitivity limit of Nqiv ~ 10^^ cm~^, which is difficult at present. They claim 
to have found associated C IV lines in 90% of clouds with Nhi > 1.6 x 10^^ cm~^ and 75% of 
clouds with Nfff > 3.0 x 10^^ cm^^. If confirmed, such a large fraction of enriched clouds will be 
difficult to explain with the two population model. As emphasized by Songaila and Cowie (1996), 
the alternative hypothesis of an enriched IGM will then be more appropriate. 

We briefly note here that a problem also prevails at low redshifts about the interpretation 
of the fraction of Lya lines with associated metal lines. Although the association of some of the 
low column density Lya lines with galaxies seems sure, it is still uncertain as to how physical the 
association is (see Le Brun et al. 1996 for a discussion). It is possible that instead of being directly 
associated with galactic halos (as in, for example. Mo and Morris 1994), the Lya lines trace some 
structure related to formation of galaxies (Morris and van den Bergh 1994, Petitjean et al. 1995). 

We note here the possibility, as shown in several recent numerical simulations (Miralda-Escude 
et al. 1995, Hernquist et al. 1995) that very low column density Lya lines may arise in mini 
pancakes and not minihalos. In this case, the density profile will be much different and our results 
may not hold. However, we note that mini pancakes can account for Lya lines with Nfji ~ 10^^ 
cm^^ (Miralda-Escude et al. 1995) and lines with higher H I column density come from the 
intersecting points of pancakes, and correspond to minihalos or galactic halos of our model. 
Therefore, our results are most certain to hold for Nhj ^ 10^^ cm~^. In any case, as explained 
above, the fraction of Lya lines with associated metal lines for Nhi ~ 10^^ cm^^ will be difficult 
to determine in the near future, and our model will not be testable for such values of H I column 
densities. 
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Incidentally, the case for an enriched IGM was described by Silk et al. (1987) much before 
the discovery of associated metal lines in Lya systems. They discussed the possibility of enriching 
the IGM by galactic winds from dwarf galaxies. The detail of such a scenario, however, remains 
to be worked out. Pregalactic objects such as Population III stars might also enrich the IGM at 
high redshifts, but again no detail models exist. 

In the face of such interesting but still poorly understood scenarios, it therefore seems 
reasonable to pursue a modified version of the minihalo model, such as the two population model 
presented here. We have shown that a fraction of ^ 0.5 of Lya forest lines with associated metal 
lines can be understood in the framework of hierarchical structure formation model, assuming that 
halos with velocity dispersion 15 ^Vc ^ 55 km/s do not produce stars and heavy elements. The 
model predicts that the fraction of clouds with associated metal lines increases slowly with the 
H I column density. We have shown that the results are robust and do not depend strongly on the 
assumptions of the limiting values of velocity dispersion. We have further used the age-metallicity 
relation of Galactic halo stars to predict the evolution of the average metallicity of the associated 
metal line systems. In particular, we find that the average metallicity at z ~ 2 should be about a 
hundredth of solar. 

We thank T. Padmanabhan, R. Srianand, and E. J. Wampler for stimulating discussions 
and helpful comments. We are also grateful to the referee for his constructive suggestions for 
the improvement of this paper. This work has been supported in part by the Grand-in-Aid for 
Scientific Research (08640318) of the Ministry of Education, Science, and Culture in Japan. 



A. Nhi vs. 6 in galactic halos 

We describe below how to obtain the H I column density as a function of the impact parameter 
b for Lya-absorbing gas clouds in galactic halos. 

The probability Pdl that a sightline in the range oi I ^ I + dl encounters the clouds confined 
in the ambient hot gas is given as Pdl = Ugadl, where n^ is the number density of clouds and 
o" = ttR'^i is the cross section of each cloud. Under the assumption that each sightline encounters 
one cloud in average, P is normalized as / Pdl = 1. Then, if each cloud has its own column density 
^Hiii") as a function of the galactocentric distance r, the average column density < Nhi > (b) at 
the impact parameter b may be written as, 

< Nhi >{b) = J NHi{r)nsadl 

To obtain the possible expression for Ug, we assume that the whole area of a galactic halo 
projected onto the sky is in average covered by ensemble of clouds without overlapping. This is 
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equivalent to / Pdl = 1 as described above. Then, the average number density of clouds within 
b, denoted as n^, may be approximately given as n^ ~ l/(47ri?^;6/3), and this yields the number 
density Ug at each 6 as n^ ~ l/(27ri?^^6). Equation (Al) is then reduced to 



<N 



HI 



1 rCxD TflT 



Here we note that N}{i{r) oc r^^^/^ = Kr^^^'^ aX r ^ r^ where H I column densities of Lya 
absorption lines are lower than 10^^ cm^^ in galactic halos. Substituting this into eq.(A2), we 
obtain 

Jh \/r^-b^ Jo cosh'/'^x 

where we introduce the variable x by r = 6 cosh x in the second equation. The integration in 
eq.(A3) yields 0.91, and this can be regarded as ~ 1 in view of the other uncertainties. This allows 
us to adopt the approximate equality as < Nhj > (6) ~ NHi{r -^ b) . 
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Figure Captions 

Fig. 1. — (a) The impact parameter b (Mpc) of a sightline as a function of the total mass of a halo 
M (Mq) to produce an H I column density Nhi at z = 0.5 {F = 0.05, J-21 = I, Xv = 15, and 
h = 1). The solid lines correspond to q = 2 (or J_2i — 0.18 in eq.9) for Nhi = 10^^, 10^^, and 
10^^ cm^^ from left to right, respectively. The dashed line corresponds to q = (or J-21 = 1) for 
Nhi = 10^^ cm^^. These lines are quite insensitive to the adopted cosmological parameters of Qq 
and Aq. For comparison, the observed value of ~ 160 h~^ kpc for the size of gaseous halos at z ^ 1 
is plotted as a dotted line (LBTW). (b) The H I column density Nhi (cm~^) against the impact 
parameter b (Mpc). The solid lines correspond to q = 2 for M = lO^'^, 10^^, and 10^^ Mq from left 
to right, respectively, whereas the dashed line to a = for M = 10^^ Mq. The other parameters 
are the same as those in (a). The dotted line is the observed size of gaseous halo ~ 160 h^^ kpc, 
whereas the bold solid line denotes the reported correlation between Nhi and b (LBTW). (c) The 
impact parameter b (kpc) of a sightline as a function of redshift z. 

Fig. 2. — (a) The redshift evolution of the number of Lya absorption lines with a column density 
larger than Nhi = W^'^ cm~^ (F = 0.05, J-21 = ^, Xy = 15). The range of circular velocities Vc 
are 15 < T4 < 250 km s~^ (i.e. Vi = 15 km s~^, V3 = 250 km s~^). Thick and thin lines are for 
the EdS universe {Qq = 1, Xq = 0, h = 0.5) and the non-zero Aq universe (JIq = 0.4, Aq = 0.6, 
h = 0.65), respectively, for CDM models with bias parameter bg = 1 (solid lines) and bg = 1.3 
(dotted). The different curves at z < 2.5 corresponds to the different values of a in eq.(6): the 
upper curves for a = 2 and the lower for a = 0. Filled squares denote the observed number of Lya 
lines at z ~ 3 (Petijean et al. 1993), z ~ 1.5 (Lu et al. 1991), and z ~ (Bahcall et al. 1993). (b) 
The dependence of dN/dz on the H I column density limit for z = 2.5 and bg = 1. Filled squares 
denote the results of Ranch et al. (1992) for 2 ^ z ^ 3. The other parameters are the same as 
those in (a). 

Fig. 3. — (a) The redshift evolution of the fraction of Lya absorption lines associated with galactic 
halos in the range of 55 < 14 < 250 km s~^ (while 15 < 14 < 55 km s~^ for minihalos), with 
Nhi > 10^^ cm~^. The error bar denotes the Keck results for the fraction of Lya lines associated 
with CIV metal lines (Cowie et al. 1995; Tytler et al. 1995; Womble et al. 1996; Songaila &: 
Cowie 1996). The meanings of the four curves presented are the same as those in Fig. 2(a). (b) 
The dependence of the fraction of Lya absorption lines associated with galactic halos on the H I 
column density limit at z = 2.5. The others are the same for (a). 

Fig. 4. — The dependence of the fraction of Lya absorption lines associated with galactic halos 
on the lower bound of the circular velocity V2 (while Vi = 15 km s~^ and V3 = 250 km s"-*^), for 
z = 2.5 and Nhi > 10^^ cm~^. The meanings of the four curves are the same as those in Fig. 2(a). 
The observed fraction of Lya lines associated with CIV metal lines, 0.5 ~ 0.75, is indicated as a 
bar. 

Fig. 5. — (a) The redshift evolution of the average metallicity < Z > of absorption lines (solid 
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curves for 6^ = 1 and dotted for bg = 1.3), based on the model described in subsection 3.3. The 
range of circular velocities for galactic halos are 55 < Vc < 250 km s~^. For reference, the metallicity 
evolutions if the star formation commences from the beginning of the universe, tsF = 0, are shown 
as dashed lines. The crossed error bars denote the Keck results for the Lya lines associated with 
CIV metal hues ([C/H]) (Cowie et al. 1995; Tytler et al. 1995; Womble et al. 1996), and filled 
circles denote the abundances of damped Lya systems ([Fe/H]) observed with the Keck (Lu, Sargent 
& Barlow 1996). (b) The dependence of the average metallicity of absorption lines on the lower 
bound of the circular velocity V2 (while V3 = 250 km s"-"^), at z = 2.5. 



23 



(a) 



-1 



o 



O 

^ -2 













- 






//// 


- 


- 






//// 


--- 








■y/ -y -V y- 








/ / // 




- 


NHr 


-XQi^^cmryX 


/// 


- 


- 






/ 


- 


- 


X^. 


X^ 


z = 0.5 


- 




.Ay 


Y , , 


a = 2 

1 1 1 1 1 1 


- 



6 

o 



o 



16 



15.5 



15 



14.5 



14 



10 
Log M (MJ 



1.5 -1 

Log b (Mpc) 



12 



-0.5 



14 



(b) 
















\ z = 0.5 
\ a = 2 


- 


\ 


V 


\ \ ^ 


\ \ 


- 


- 


\ M=101°Mq 


\lO"\ 


\ \l012 


- 


- 


\ 


\ \ 


\ \ 


- 


- 


\ 


\ \ 


\ \ 


- 


- 


\ 


\ \ 


\ \ 


- 


- 


\ 


\\ 


\ \ 


- 




1 1 1 l\ 1 


1 1 1 1 \ A 1 1 1 \i 1 \ 


- 



500 



24 



400 






a = 2 



M=10i2M, 



o 



300 



200 



10>'M, 



100 



lOi^M, 



lOiiM, 



0.2 



0.4 



0.6 



0.8 



25 



1000 



S 100 =- 



13 



•z. 

X! 



10 37 



•z. 

A 
N 

\ 



1000 



100 =- 



Nj^,>10i4 cm-2 




14 16 

Log Nh, (cm-2) 



26 



(a) 



0.8 



0.6 



0.4 



0.2 



Mini: 15<V^<55 km s^' 

Galactic: 55<V<250 km s"' 



Nh,>101* cm-2 




(Ho ,Ao ,h ) 
(1 ,0 ,0.5 
(0.4,0.6,0.65) 



(b) 



10 



1 + z 



Mini: 15<V^<55 km s-' 

Galactic: 55<V <250 km s-^ 



b =1 



z = 2.5 




0.4 



0.2 



(Ho ,Xo ,h ) 

■(1 ,0 ,0.6 ) 

(0.4,0.6,0.65) 



_L 



12 



14 16 

Log Nh, (cm-2) 



27 




120 



V (km s-i) 



28 



(a) 



Njj,>10i* 



tgj, from Iff 



(Ho Ao ,h ) 
(1 ,0 ,0.5 ) 
(0.4,0.6.0.65) 




1 + z 





ISI 



A 
CO 
V 

M 
O 



(b) 



- 



1 - 



I ' ' ^^r^^ ^^ 

(fio ,\ .h ) 
^^^ (1 ,0 ,0.5 ) 

^^^ (0.4,0.6,0.65) 




40 60 80 

Vg (km s-i) 



100 



120 



